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Boday,  deterministic-differential-equation nmodeis( that 
are commonly called Lanchester-type combat models) are 
widely used in military. operations research for defense 
planning purposes, especially as concerns land combat. In 
conbat analysis, such questions frequently arise as ‘Who 
will win the battle ?, How long wiil tne battle last ?, How 
Many survivors will the winner have ?, and How do the force 
levels change over time ?',etc. There have been a lot of 
studies to estimate attrition in battle to answer such ques- 
tions as these. 

In 1914,F.N. lanchester (Ref. 1] first proposed a mathe- 
matical model of combat in order to justiiy the princirle oZ 
Bomeentraticn of forces under "modern" conditions. de 
Lelieved that in "ancient" times war was a series of 
one-to-one duels between nen engaged in hand-to-hand conkat 
Hsing sword, axe, larce, etc. So tnere would be no advan- 
tage gained from concentration of forces. Lanchester 
reasoned that modern technology nad changed the nature of 
this "ancient" warfare and now there was a decided advantage 
meme concentration of forces".1 He hypothized that each 
side's casualty rate was proportional to the number of 
ESposing firers. As LEA Out noc. yer explicit) in 
EcHester's original work is the concept of an "attrition- 
Fate coefficient" as the rate at which a single firer kills 


a particular enemy target type. Furtsermore, Lanchester 


_1 Helatold notes that the analysis of the. data from 92 
NE courat battles_ sugyests_ that victory in battle is 
Exmarily determined by the factors other than numerical 
ENDETTOLI1ity and challenges the abiilty, os: any model of 
compat which concentrates almost exclusively on numerical 
ENSSC size to yield a practically useful predictor of 
Story in battle. 





Eught that the number of survivors were the critical 
ED or deciding the outcome of the battle in "modern" 
warfare. 

Since Lanchester's original efforts to describe the 
dynamics of combat mathematically, many analysts - in an 
metempt to add reality to the combat description - have 
extended the theory to include additional factors. Many o£ 
these extentions were described by Dolansky [Ref. 2]. In 
his paper, Dolansky pointed out that the  Lanchester 
attrition-rate coefficients were aard to determine for 
particular weapons, and, accordingly, coefficients of known 
Magnitude had been assumed in most aodels. Thus,  Dolansxy 
concluded that the use or  Lanchester-type models for 
prediction or battle results had been hampered by this 
Mon lity to predict rumerical values for attrition-rate 
coefficients. However, even if we nay pe in doubt as to the 
proper rorm of attrition relations, we may at least te aple 
emake rragmatic assertions such as, 'If tke Lanchester 
attrition 1s used for predictions, then by measuring certain 
MAS tors and performing certain operations a prediction ol 
AA anba=such an aspect of combat can be obtained.' So, 
EBthouch there has been a continuing discussion among mili- 
tary operations analysts abcut the merits of Lanchester-type 
attrition-rate, it is still generally accepted as the tasic 
founding for all the combat models studies and attrition- 
rate coefficient is still a core of these studies. 

Bonder/Sarroot { Ref. 3] proposed that an attrition-rate 
coefficient be defined as the reciprocal of tne expected 
time for anindividual firer to kill an enemy target. A 
more thorough dicussion, however, Oe Eno Sti fication For 
Eonger/Sarefoot is to be found in Taylor (Rei. 3]. Barfoot 
(Ref. 5] suggested that the attrition process be considered 
Alsa. rereval process. Pus ;esOs a Singke Y firer, we will 


Bine that each individual Y kills X targets according to 
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an attrition process in which times to kills are indepen- 
Mentiy and identically distributed (1.1.d.) random variables 
E.V). We will also assume that each firer acts indepen- 
Bently OL any other firer, and Sinilarly for targets. Suen 
an attrition process turns out to bea renewal process in 
fee theory of stochastic processes. 

Once we introduce tne theory of renewal process, we can 
consider many other random variables such aS tae time up to 
en kill, S; the number of kills in time t, N(t); the 
expected number o eS In tiren, A(t), etc. Thus this 
thesis starts with consideration of Bonder/Taylor's concept 
of attrition-rate coefficient, aad then studies the prob- 
du lity distribution of the tine between casualties 
according to the thecry of renewal process. 

In this thesis, we will mainly focus on the fundamental 
Lanchester-type attritıon paradıga waich deals with homoge- 
neous forces, because understanding this basic paradigm is 
essential fcr extentions to further coaplex models. 

HUS thesis is organized in the tollowinc fashion- 
First, we review a few relevant facts about the Lanchester 
mereition modeis. Then we review basic ideas of renewal 
process and develop the random variables whose study is of 
our concern, followed by the study of these random vari- 
aples. Then we will describe tne combat attrition process by 
Semparing the values of different random variables obtained 
througn renewal process. Finally, tae conclusion mainly 
contributes to the summary of tne work and suggests further 


Studies. 
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II. FUNDAMENTAL LANCHESTER-TYPE ATTRITION PARADIGY 


Àe THE BASIC PARADIGM 


Let us consider combat between two homogeneous forces: a 
homogeneous X force opposed by a homogeneous Y force. We 
will focus on the force-on-force attrition process in tae 
compat ketween these two homogeneous forces. The  rasic 
Lanchester-type paradigm for "modern" warfare assumes that 
the casualty rate of such a homogeneous force is directly 
proportional to the number oi enemy rorces, e.g. the X force 


casualty rate is given by 


dx ou _ 

Fs =i= ay dith X(Q9) = > $ 
and the Y force casualty rate is given by 

D CM aith yO =y , 

dt 0 


where ‘a’ denotes the rate at which a single Y rirer kills 


A ER CHER cub  UNEP SS > O: ls CERO O. A mm m 


cient. Here Se and UN denote tae number o£ X and Y 
compatants (respectively) AE ne De jinnin or the catcle. 
Many factors may affect tke value of 'a' and 'b', but for 


the time being it is not essential that we be explicit about 
functional dependence or a and b. Let us next address the 
problem of computing a reliable numerical vaiues for such 
attrition-rate ccefficients. 
Two approaches have been developed in the United States 
to determine the value of these coefricients as follows: 
(1) à Statistical estimate based on "combat" data 


generated by a detailed Monte Carlo combat simulation. 
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Cana ttcal surmodel of attrition process for tae 
Bor ae OOND a naton or $orer and target type. 
In this thesis we will consider oniy the second approach. 
ENGonder [ Ref. 6] has called this approach the use © 
free-standing or independent analytical model. The kasic 
idea of this approach is to develop an analytical expression 
for each required kill rate by considering tae single firer 
engaging a 'passive' target (i.e. ome that doesn't fire 
back) and then tie all the attrition rates together in 
force-on-force combat with a Lanchester-tyoe model. 
px .ootolTRet. 5] in his study bas suggested taking a 
Lanchester attrition-rate coefficient as the reciprocal of 
we expected time for an individual firer to kill an eneny 
“target.'? Within ccntext of the above homogeneous force 


Lanchester-type coanbat model,this means 


Ec; denotes a r.v. representing the tiae for a single 
Y EE. to kill an X target and E[T] denotes expected time 
to Kili. Taylor also supported bSbarfoot's idea in his recent 
study [Ref. 7] by ccmparing times between xills that were 
both exponentially distributed and non-exporéentially 
distributed, and using previous equation to determine the 
expected time to kill a target. 


“However, , the justification of this approach is not 
accepted and is No E sou somewhat controversial. But in 
Case of exponentially distributed times between kills, a 
deterministic Lanches cr E model may indeed be considered 
Nela tae mean course ol Combat. Anyway, this approach is 


ESE pted as generai guide for further study of attrition- 
Bare coefficient. 
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A DETERMINATION OF EXPECTED TIME TO KILL 


Bonder aná Farrel (Ref. 8] have developed general neta- 
ERN cro derterun?n:nLg the attritior-=rate coefficient by 
determining the expected time to kill a target, z[T], fora 
wide spectrum of weapon system types. Morever, research 
since the mid-1960's has led to the developzent of several 
other methods for computing the expected time to kill a 
target. For present purposes, two methods will be of our 
interest, 

(1) Method based on sum of component event times, 
and 
(2) Method based on first-passage time in Semi-YJarkov 
process 
If we assume tne special case of tactical interest: namely, 
the case of Markov-derendent  fire,Bonder [Ref. 3] has shown 


BOSE C R[T] turns out to be as rollows: 


(t +t) (t +t) 

E[T]et * t,-t a a e ia + P(H|B) - P) 

a 1 un P(K]1H) P(H[ 4) P(KI&) 1 
where all symbols are defined in Table I This expession for 
Mee! holds for the fcllowing assuaptions; 

(1) Markov-dependent fire with parameter P , P(H]H), 
and P (H|) 

(2) Geonetric distrioution for tie nunber of hits 
required tor a kill with parameter P(K{d) 

En assessing the tue ractor to determine tae 
Eu tlon-rate coeificient, it turned out that target acqui- 
sition process and weapon system capaoility were most impor- 
tant factor from the foraula. However, when we consider the 
time factcr of the attrition process, tame acquisition 
process would be the most important factor to study because 
ENUrectly affects the required lengtk of time to killa 
target. So we will elaborate on target acquisition process 
for further use in attrition analysis. 
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TABLETI 
Variables for Expected Time to Kill 
time to acquire a target, t 
time to fire first round after target acquired, t 
time to fire a round following a n. 
mime to fire a round following a niss,t 


m 


meme Or flight of the projectile, 7 


ti, 


Esobatility ot a Eit on first round, E. 
Peobatitity Of a hit on a round following a hit, P(HIH) 


Mmeobabitity Of a Lit on a round following a miss, P(H] 1) 


EpEobabllity orf destroying a target given hit, P(X[H) 


A T———————— —————————M——— Lo a rS | mannan! 
. 


| 


TARGET ACQUISITICN PROCESS 


Target acquisition is the initial step in engaging 4 
target and a most deciding factor that affects the time to 
fees AN important distinction made in Vki'ts Lanchester-type 
combat models is whether the target acquisition process ol a 
ENngle "typical" firer type isa serial process or a 
parallel process [Ref. 4]. The two modes for the target 
acquisition process considered by VRI's models( including 
METOR-2 ) are as follows: 


(1) Serial Acquisition 


Mey Parallel Acquisition 
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Here a single firer using serial acquisition does not 
acquire targets while engaging another target, say target A. 
When such a firer ceases to engage a target A ( due to kill, 


lost or any other reason), then ne must acquire a rew target 
all the way fron the beginning. So it is assumed that he 
does not remember any acquisitions made prior to engaging 
the target A. On the other hand, a firer using parallel 
Mecuisition searchs for targets continuously, "even while 
engaging target A and remembers those targets that have been 
acquired. Thus when such a firer finishes the engagement 
with target A, he can inmediately saift engagement to the 
next target, provided that such a target was acquired and 
processed during or before the engagement of target A. For 
both methods, the VRI models assume tnat a firer never 


engages a killed target again. 


MEN serial Acgquisition 


Now when we consider the basic Lanchester-type 
combat nodel in the previous section, a firer in serial 
ÆC ess mode or target acquisition starts again ¿ron the 
beginning in searching a new target wherever the previous 
target nas teen killed. Thus we assume that tarcet acquisi- 
tion tine is always needes to accuire a diizerent target for 
engagement. Here itis assumed that the total-force kill 
MWE 1S just the single-Zirer-xill-rate times the nuaber of 


Brers, e.g. 


dx 
dt 


E = ay, 
since each firer is independent from any other firer. Aliso, 
in the case of a serial process for each firer, the expected 
time to killa target is the same as has been shown in the 
equaticn for E(T] in section B. A Summary of attrition-rate 
coefficient results for serial acjuisition is given in 
Ile II. 
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TABLE 11 


Attrition Process Using Serial Acquisition ode 


IS 


E[T 
E 


A IO a Y firer to kill an X target 


-anm a A a e D E a D eee Te TG Te ze 


(eet otatte) 

h f a IA) 

E[T] = t #t -t t--—----t--—--—---(---------- 4 
a 


?P(H|H)-P ) 
1 h P(XIH) P(HIM) P (KIH) 1 


CN 
E 


ho oros —— —— 


As we  menticned earlier, a firer in the parallel 
mode of acquisition continues to acguire new targets even 
Mide he is engaging a particular target. Once this target 
has been killed, then he car immediately shi£zt rire to a new 
target which has been acquired aiready. In this case, 
eral- force kill rate is given bj the product of the kill 
rate of single firer against acquired targets and the 
expected number of  firers who have already acquired one or 
more targets. 

Ihe summary results for attrition process using 
parallel acquisition is given in table III . Note here that 


E is same as the vrcbability that a single Y firer using 
XY 


1 





parallel acquisition mode has available one or more acquired 


meegets at which to fire at time t, because if ne is firing 
Et 


we know that he already acquired the targets. Then 


should 


also be noted that 1 y represents the expected number ož Y 
XY 


firers who have already acquired one or more X targets. 


MW ———— 


Attrition Process Using Parallel Acquisition Mode 


R : 
xy 


TABLE III 


Propisingle Y firer using »arallel acquisition 
S ceing at A target at random tine t } 


ooo — ove sek (s)asj 
-0 xy 


IDE: kill pate against acquired targe 


IN OE aire to KI) an acquired X target 
COLL i tional kill time) 


1E WE eea aca r Cirer aeguires X targets 
Ie t Ver appropriate 


CE tt 
Eu un ) 


= t -t + ------- + ------- (7--------- *P(H|H)-P.] 
1 h PI P(HI1)  P(K|H) 1 
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Now we may assume that it would take longer to xiil 
target for a £irer using serial acquisition than a rirer 
using parallel acquisition because of the time fraction of 
BE anisition of target. This distinction between tnese two 
acquisition modes may be further studied as the case ofa 
delayed renewal process. Here we assume that expected time 
to kill the first target wouid be optained from the serial 
Mmegrisition but following tine to kill would be fron the 
parallel acquisition, provided that the targets are acquired 
already. But in this thesis, we will only deal with the 
firer who uses the same acquisition mode to attrit all 
purgets. 

Before continuing, it is important that the reader 
understand both the basic  Lanchester-type attrition model 
and also Bonder's approach of the attrition-rate coefíicient 
as the recicrocal of expected time to kill a target. It nust 
te recognized that this expected tine to kill a target is 


eerimarily dependent cn the accuisition process. 
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Tit. ANALYSIS OF THE COMBAT ATTRITION 


A. INTRODUCTION 


The Lanchester attrition-rate coefficient is the rate at 
which a single firer kills a particular enemy target type in 
Lanchester-type conbat. Development of technically-sound 
and scientifically-valid methodology for determining numer- 
ical values for Lanchester  attrition-rate coefficient is an 
essential prerequisite for building a Lanchester-type combat 
model. The tasic construct cr B3onder/3arfoot methodology was 
to take a Lanchester attrition-rate coefficient as the 
Mecitrocal cf the expected time for an incividual firer to 
kill an enemy target as mentioned in the previous charter. 
May lor provided justification (Ref. 4] for taking Lancnester 
attrition-rate coefizicient as the reciprocal of the expected 
mme to kill He started tren the basic hypothesis that 
combat is a complex random process, but it contains enough 
regularity that the appropriate Lanchester-type equations 
meena gOCd approximation to the mean course of combat. It is 
clear that the casuaity rate 1s equal to the reciprocal of 
Mier expected time for a force to inflict a casualty, when 
the tines between casualties are exponentially distributed. 
However, in case the times between casualties are no longer 
exponentially distributed, Taylor used the suggestion nade 
PyesOnder and Barfoot. Bonder [Ref. 3] and Barfoot [Ref. 5] 
Suggested defining the Lanchester attrition-rate coefricient 
as the expected rate which a single firer kilis enemy 
targets. 

the Spirit of Bonder and Barroot, Taylor provided 
Ec rigorous  justizication (Ref. 7] for the Lanchester 


HM Utlon-rate coefficient that does not assume an 
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Pe ogential distribution Eor times between Casualties. He 
considered the case where the initial force size of xX and Y 
EN .- enousa t9 insure a negligisie probability that the 
battle is terminated before an attrition pattern coula be 
established. He made no specific assumptions about the 
distribution oí times between kills, but assuned that each 
My dual Y force kills X targets according to an attrition 
process in which the tines between kills are independently 
and identically distributed random variables. Thus in the 
parlance of the theory of stochastic processes, he said that 
mich an attrition process is called a renewal process. 
Prior to this, attrition prediction hnas been difficult due 
memtne inability to ¿redict casualty patterns. By investi- 
gating the distribution of the times between casualties via 
renewal thecry, combat attrition analysis becomes easier. 
Indeed, we know a great deal anout tne casualty pattern once 
we know the distribution of the interarrival time of kills. 
AS we consider more variables fron the renewal process, we 
ind aore information about tae casualty pattern. 
Furthermore, we can track the development of the contat 
Beer tion more precisely by studying the probabilistic 
distribution function of those random variables that can be 
considered in the renewal process. 

we assume here that the reader is aware of the concept 
Of a counting process. nen we think about combat between X 
and Y forces, we can assume that tine between kills bya 
single firer has sone distribution F. Now we are interested 
in observing the occurence of casualty and the number of 
Asta lties, N(t), that have nappened in time interval (0,t]. 
Furthermore, when times between Kills are i.i.d., we call 
this counting process as renewal process. 

In this chapter, NS INEO attrition atalysis will 
include the foilowing : 





Me Ev.) Nit), total number oc casualties by tine t, 


it's distribution, and properties. 
e The expected number of casualties by time t. 
e Theorems discussing casualty occurence. 
e The limiting behavior of the casualty occurence. 
e The age and excessive life or the casualty occurence. 
e Examples with specific casualty distributions. 
e Application of renewal process for casualty estimation. 


total force attrition by pooling the single firer 


renewal processes. 


Be REVIEW OF RENEWAL PROCESS 


Ii the sequence of nonnegative randcem variables (X ,X , 
EN are i.1.d., then the counting process (N(t), t20 j is 
said to bea renewal process (or oftea called an ordinary 


renewal process). Thus, a renewal process is a counting 
process such that the tine until the first event has scme 
distribution F, the time between the first and the second 
event has, independently of the tine of the first event, 
the same distribution F, and so on. So when a event occurs, 
we say that a renewal has taken place. 

An example of a renewal process, let us suppose tnat we 
have an infinite supply of lightbulbs whose lifetines are 
l.i.d.. Suppose also that we use a single lightouib ata 
time and when it fails we immediately replace it with a new 
one. Under these conditions, {N(t), t20 } is a renewal 
process when N(t) represents the number of lightbulbs that 
have failed by time t. Accordingly, if we assume that tnis 
peere is a kili by a single firer - X or Y -, then killing 


process also can be considered as a renewal process. Thus, 


2? 





renewal (or event or arrival) will be synonymous wita kill 
Dr casualty occurence in the rest of this thesis. Je now 
assume that casualty cccurence is a kind of renewal process. 
me, in particular, the distribution of interarrival time of 
ESI is exponential with probability distribution function 

E W then this renewal process turns out to be a Poisson 
process with rate, . [Ref. 9] 

Thus we may assume that X casualty occurence due toa 
single Y firer who uses the same acquisition node to engage 
each X target can be explained by ordinary renewal process 
because we assume that times between casualties are ail 
i.i.d. However, we often consider a counting process for 
which the first interarrival time has a different distribu- 
tion than the remaining ones. So, we may think of a Y firer 
who at first uses serial acjuisition mode to engage an X 
em cet, but once he has finished engaging an X target, then 
he nay use parallel mode to engage tne rest or tne targets. 
So, in this case, we can assume taat interarrivai tiae of X 
casualty has different distrikution between the first and 
the rest. 

Foraally, let [X ,i21) be a sequence of independent r.v. 
witb A Eom det Dion Aman 22 naving distribu- 
Mon F which is different [rom G. Then counting process 
[N (t),t2) is said to ce a delayed renewal process. SO, we 
can think of using delayed renewal process for different 
acquisiticn mode, but here in this thesis, we will only deal 
with the case of ordinary renewal process for sinplicity. 

we now define for the renewal process a number of asso- 
ciated random variables whose study is the objective of this 
thesis. Table IV and Figure 3.1 give the intuitive interpre- 


tation of the random variables obtained in renewal process. 
Here we can interpret each random variable as roilowings: 


e X is the time between the (n-1)st ani n-th X force kill 
n 
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TABLE IV 


Random Variables in Renewal Process 


called as spread 


| 
| | 
| E. V interpretation | 
| VE 
time between the (n-1)st and n-tà renewal | 
| a i-e. n-th interarrival tine | 
| R tine up to n-th renewal | 
| - | 
N (t) total number or renewals in (0, t] | 
| M (t) expected number of renewals in (0, t] | 
| known as renewal function 
eccL E e | 
Zt length oí time measured fron the last | 
renewal until a given tine t , calied as | 
age or backward recurrence time | 
p uM. . NM C ee mwen | 
YE) length or time measured fron tine | 
| t until the next renewal, called as 
EEUU eee oe cae tine 
| Sit) length of renewal time at tine t, 
| 
| 


| tame nana. Dro caa im amma ie ge, Oe nz m ET eee tocan 


B oo 


kwa single Y firer. 
moe aS the time up to n-th kill. 
n 


ENU(t) 2s tne total number of kills in (0,¢}. 


24 





i 


cas 


IA 572 (8) --> IST) --> 


a nr a a en ae a ca 

A 

i 

i 

l 

i 

l 

l 

l 
un 
pt 
ct 
eo 

| 

i 

i 

i 

i 

i 

i 
a 

Der at, Tes a. octet Meee que ras, A Sete a BS oh ra nm 


| time 
ee Tro er mn pen SS... e. e X EE cmm | TM S 
| S S S E S 

1 2 N (t) N (t) *1 
0 

Figure 3.1 Benewal Process. 


e M(t) is average number oz kilis in (0,t]. 
e Z(t) is the time between the last kill and a given time t. 
e Y(t) is the length of time from time t up to the 
next kill. 
EE (t) is tae length of lizetime of X force soldier at 


tine t. 


in this thesis, we will oniy deal with X force casualty 


EN firer to avoid the confusion ia notation. 


C. DISTRIBUTION OF TOTAL NUMBER OF CASUALTIES 


Suppose tnat for az ordinary renewal process a single Y 
firer uses the parallel acquisition aode to engage passive X 
targets trom tne beginning of combat, when we assume that 
targets are already acquired. Pret rie. CiStr pation or 
N(t), the total nuaber of X force caSuaities by a single Y 
Beurer, can be obtained at least in theory by first noting 
the relationship that the number of kills by time t is 


greater than or egual to n ii: ard oniy if the arth kill 


Occurs before or at tine t. 


Formally 
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Men «se» S «t (3.1) 
n 


Fon 3.1, e octain 


P(N(t)=n) = P(N(t)2n) - P(N(t) >2n+ 1) 
= P{S <t} - P{S st) 
n n+ 1 
Sew, since the random variabie X "s, (i21), are independent 
1 n 
aná have a common distribution F, it follows that S = =. 
DNE 


=1 1 

MW -tributed as F , which is n-zold coavolution oz F with 
n 

itself (Ref. 10]. Bere eee tne Probability distributio 


e tine between casualties. Thererore we obtain 


P{N(t)=n} = F (t) - F Ca S t) = 1 
n n+ 1 0 
and 
UNE) En) = P(S 2t) = 1 = Fr (t 
w(t) Sn} { a j M ) (3.2) 
BEnererore the probability distribution of i(t), which is 
total number oi casualties by sirgle Y firer, can be 


ErGuned explicitly fcr all n. 

The simplest special case of ejuation 3.2 is obtained bv 
taking the renewal prccess to be a Poisson process. 
Then S has the special Erlang distribution with n stages 
because each Y is exponentially distributed. Also, it is 
known from the previous argument that N(t) has a Poisson 


distribution with mean t. 


Bee EXPECTED SUMBER OF CASUALTIES 


des Renewal Function 


nD A O CONES ae 2 a ES SS = 


Expected number of casualties, E[N(t) ], which is 
defined to be the mean number of casualties in the interval 
(0,t], 1S known as the mean value or renewal function, M(t). 
From equation 3.1, the expected ruaber of casualties in tine 


— i 


Bean be obtained as following: 
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NAS 


coe. 


2t) AE) J = 1 


K 


D) 


=> 2. P{N(t) =k} 
n= 1k=n 


MY 


ta 


=> _P(N(t)2n) 
n= 1 


bo 

NES 
n= 1 n 

haz 
PEA E) (3. 3) 
n=1 qn 


But here it is generally very difficult to find M(t) 
Epectly from the distribution F. Thus, we will employ the 
Biblace Stieltzes Transform (L.S.T.) (Ref. 9] to find more 


conveniently an expression for M(t). we generally need 
three steps to find M(t). 


Step 1: 
Fror tie definition orf L.5.T,., we know that L,S.T. 


e (t) is 
Dy 


* | -st 
M(s) un c dM (t) 


and also 
* ea E 
F(s) =i e SEES 
Step 2: 
Ihen from equation 3.3 
o9 
* - 
Mis) =~ F (s) 
n= n 
(“a x n 
=> “P(s)] for s»0 
n= 1 


The alove is a sum of geometric series and can be written as 


follows: 
ez 
on 
1-F(s) (3. 4) 


D 





or equivalently 


x 
* 4 (s) 
F(s) = uos 


1*M(S) 


tep 3: 

Since the L.S.T. of a function uniquely deternines 
the function, the renewal function deterzines the distriru- 
tion of the kill time and therefore determines tke prob- 
ability law or the casualty occurence. Specifically there is 
a one-to-one correspondence between interarrival tine 
AS tribúution F, and renewal function, M(t). So if we have 
the distribution F, then we can find A(s) from the equation 
EN. But to find 4(t) we nave to transrorm M(s) to H(t) by 
ENWerse application of L.S. T. 

Let us try the ezanple. Remenber that this exanple 


will be applied all the way through this chapter. 


Example 3.1 


ENbpose that sach interarrival time to kill X,,(i»1) nas 
i 
EEorLrential distribution with . Inen L.S.T. of F is 


x 
+ t -st_ 
ts) = e aF(t) 
20 
CER Sc ESSE 
= | e ^e dt 
- 0 
F -(2+s)t 
AA, - (Ars)t = 
Free (^*s) 
AtS 0 
Dx 
= ----- if s>- A 
ARS 


Men £rom equation 3.4 
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IEA ET 
* A-(VIA+S)] S 
1=F(S) 


where 1/S was an L.S.T. of F(t), which had value of t.3 So 
M(t) =At + c (some constant), but we know that M(0) = 0, 
because there would be no casuaity at time O. So, finally 
me) = At. 

Let us be more Specific about this example. Suppose 
that interarrival times ot: casualties are i.i.d. and expo- 
nentially distributed with zean 30 minutes. Then 

a) What is the distribution or total X casualties at 5 
hours later ? 

ENS hat is the probability that 3rd X kill was occured at 
time t = Y ? 

eye What 15 the prckrability that at time t = 2, no X casu- 
EN occured ? 

To answer these questions we know tnat M(t) - 2t from above 

exampie because = 2, then number or X casualties in time t 


Bas Poissch distribution with mean 2t. 


3Let F(t) = t for t20. Thea L.S.T. of F(t) is 


* SE (Bs -st 1 -st 7? 
F(t) = ^, e ar (t) = e AE KES TI ===- e | 
PEU 0 S 2) 


1725702520 
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Tnen 


-10 n 
€ (10) 
DIS =D = ES 
nu 
3 -8 
8 e 
b) P{N (4) =3) = er = „0235 and 
0 -4 
4 e 
c) PIN (2) =0] = o = .0183 


Now an integral equation for M(t) may be obtained by 
Bonditioning on the time of the first casualty occurence. 


Upon doing so, we obtain 


YO r 
“(t)= 2 EN (t) IX > NI) (San 
however 
0 O E 
Bent) IX = x] >= i| 
1 nai + M(t-X) SE (375) 


Er: the first casualty occurs at time x, xst, then fron 
this point cn the process starts over again, and thus the 
expected number of casualties in (o,t] is just 1 plus the 
expected number of casualties in time t-x from the beginning 
Mana equivalent renewal process. Putting equation 3.6 in 
3.5 yields 


at 
M (t) ie - 1*M(t-x) JdF (x) 


FE 
ESEYE x MIU E SE) (223) 
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where equation 3.7 is known as the renewal eguation anc r is 
a known function and M is an unknown function to pe deter- 
mined as a solution to the integral equation 3.7. ROSS 
providei the solırlonssromsmeguatıon 3.7 [Ref. 11] as 


following and the solution is accepted true for every case. 


Expected number of casualties up to tiae t is 





M(t) » F(t) * l. F (t-x) dM (x) (3. 8) 
where M(x) = = F (x) 
nz1 n 


EU t allows to find M(t) using oniy the known distribution 


function F, which is the distribution of time between casualt- 


NECS thus avoiding finding individual F. 
n 


3. Renewal Density 


Now let's define the renewal density as m(t) = 
di(t)/ydt. Remenber that M(t) is the ezpected number of casu- 
alties up to time t. So a(t) nay be regarded as the expected 
number of casualties per unit time; more precisely, the 
average numkec of casuaities during tne time interval fron t 
BEP tth is 


/t+h 
M(t+h) - M(t) = 3, m(x) dx (3.9) 


Differentiating tne renewal equation, we obtain the renewal 


equation for the density nm(t) 


(t 
m(t) = f(t) $t 


rea cu xq rop t»9 (321) 
So,tais knowledge of the density £ of interarrival time of 


casualties is sufficient to determine the density m(x) of 


casualty occurence. 
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E. LIMIT BEHAVIOUR OF CASUALTY OCCURENCE 


ke know that total number of casualties, M LN go to 
infinite as time goes to infinite, when we assume that force 
Size is infinite. But it would be useful to know at what 
rate N(t) approaches infinity. That is, we would lixe to be 
able to say something about limit of N(t)/t as time goes to 
feria nite. 

From the definition or E.V. we know that 


S Seo Ss 
N (t) N (t) *1 


wnere S must be the time or the casualty arter time t. 
N (t) * 
By refer to [Ref. 10] we know that 


N (t) 1 
a ooo dots 0 (8:45 
t u 
where u = E[X j = | a MIS Eo waich is the 
l 


expected time petween casualties. 


me the function 1/u is often called as the rate of the 
renewal rrocess. Thus the average number of casualties per 
unit time converges to 1/u . Then, how about the expected 
average number o£ casualties per unit time? Is it true that 
Prey /7t also converges to 1/u as t -=> ? This result, known 
as the elementary rerewal theorez, will be stated without 


proof (Ref. 11]. 


M (t) 1 
------ --->  --- as t --»5 M (3. 12) 
E u 


Thus whatever the distribution F is, the average nunbe 





t 
M(t) = --- (3. 13) 
u 


Also, the average number of casualties in the interval 


OK u AA mmm A-—A A A ES ai am m En EG ao a GENS 


(t, t*h) for h>0 is approximated for large t as following, 


h 

M(t*h) - M(t) = --- as t --> (o (3. 14) 
u 

for any distribution F with mean u. 


Here equation (3.14) is known as plackwell's theorem when F 
is not lattice. 

From the example 3.1, where casuaity occurence time is 
exponentially distributed with mean 30 minutes, the average 
number of X casualties up to time t is approximated ry t/u, 
Mich is 2t. Then after 2 nours or combat, we can assume 
that the average number of X casualties by singie Y firer 
fet be 2t, which is 4 casualties. 

Another limiting result which may prove to be useful in 
rerewal process also concerns tne r.v. N(t). It is shown 
[Ref. 12] that N(t) has an asymptotically normal distribu- 


tion with mean t/u and variance t 7/u'.  Ihus 


Eom PIN(t) Sy) =p zen (2,905) 
DD xal 


Suppose of the same example 3.1 that casualty occurence nean 
time is 30 minutes. Then u = 0.5 = B[T] = 14 . SO = 2. 
eo Var(T}] = 1%: = 1/4 = 0.25 hf , which is 900 min: here. 
NS assume that we had a 10 hours or combat action. Then 
Eupectea nurber of X casualties is ELN(10)] » t/u = 10/0.5 = 
20 and var[N(10)]= to'yuè= 10x0.25/(0.5j= 20. 

Suppose of the question that how many X soldiers will be 
Seated in order to with probability of 95% that at least 
Mere than cna X Soldiers wili survive after 10 hours of 


compat ? This Juesticn may ke answered as following, 
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0 95 í (10 L 3 «— | A = 20 \ E goo me 20 B 
. EEEN Kae A TIT 


Mow from the normal table, 9.95 percentile of normal distri- 
bution is 1.645. So 1.645¢ (x - 20) /4.472, where the answer 
turns out to be x227.36 . So we can say that X force needs 
at least 28 soldiers as an initial force size to ke 95% sure 
that at least more than one X soldier can survive arter 10 
hours cf combat. 

There also exist a alternative normal approximation, 


when we assume that for integer n, 
P{N(t)<n} = Be where S = = X 

Brom equation 3. 1. 

Then P{S >t} SEM P {S St) = 1 - È! -------- | (3. 16) 


ODO we assume that 5 has normal distribution with ıean 


n 
nu and variance nj. II we apply equation 3.16 to previous 


exanple, 
— 110 - 0.5x \ 
0.95<P{N (10) <x} = 1 - 2 ------- ---- 
e 0.5 y x E 
En WENDE 10 - 0.5x 
Monn yo = =| Ss 0.05 or  --------- ££ - 1,645 
oss UNE 


IE SOLVE for x, it turns out to ne x228.83. So answer is 
almost same when we assume some roundiag errors. 
USAGE AND EXCESS LIFE DISTRIBUTION OF CASUALTY OCCURENCE 


As we mentioned earlier in this chapter, we can also 
consider such random variables as Y(t), which is the tine 


Brom certain time t untii next casualtyj occurence, 


34 





meen means residual lifetime of a soldier in combat fron 


Sereain time t if ke is stiil alive, and Z(t) , which is the 


time fron t since tke last casualty occurence (see Figure 
Be. 2). 

Nu m | 

N (t) N(t)*1 | 

+ EXC — m | 

| Z(t) (t) | 

| 

S(t) | 

| 

—€————————— nn — 


Figure 3.2 Age and Excess Life of Casualty Occurence. 


BIS means 


mt) - 5S - 
N (E) * 1 
Z(t) = io 
N (t) 
where Y(t) is called as the excess life O Na eS 


et) = Z(t) + Y(t) 
= 5 =- t+ t- 5 
Nit) 41 N (t) 
= S - S 
N(t) +1 N(t) 
N (t) * 1 N (t) 
EDD e: 
i=1 I 1=1 L 
= X 
N(t) +1 (3. 17) 
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Here S(t) is often called as spread and S(t) is also r.y. 


because X IS Toy: COO. 
ety 


Suppose or the Poisson process that X, is exponential, 
1 
Buen ELS(t) ] = ElY(t)] + E[2(t)] = 14 + 1l/a = 2/1 = 2E[X ]. 
v ZEN i 
So, S(t) can be assumed as the length or the lifetime of the 


X force soldier which is in combat at time t. 





1V. THE SUPERPOSITION OF CASUALIY OCCURENCES 


A. INTRODUCTION 


In the previous chapter, we have considered in detail 
the theory of casualty occurence through the ordinary 
renewal process of a single Y force soldier. We deal in thiis 
chapter more briefly with the superposition of several 
renewal processes, which is total X force casualty occurence 
behaviour. When we assume multiple Y force soldiers who 
Shoot at X force soldiers independently of each other, we 
may assume that all these Y soldiers kill X soldiers by 
indepenient casualty processes. So, we can say that there 
occur multiple independent ordinary renewal processes, 


simultaneous in time. 


En — — |. 





ES SI XY A e gonna 
| | ' © | 
| t t 1 
1 q g | 
| Process 2 | =>--- DLE '-——X---------- zen nn -| 
1 1 1 t 1 1 ! 
1 1 1 1 t 1 | 
| t t t t t t | 
MESS SR == == Ta tr ol ml o ooo o | 
| 1 1 1 1 1 1 1 1 t j 
1 1 1 1 1 1 1 t 1 | 
| pooled u ' t ' t ' ' ' 
| output IM am on aa adn ME ceres yan ey 
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Figure 4. 1 Superposition of Casualty Occurences. 
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Suppose that we have ‘pt independent casualty processes 
in operation siuuitaneously, all with the same P?.D.F. or 
casualty occurence time. Consider the sequence of casuaity 
occurences formed by pooling the individual processes. jus 


4.1 illustrates the special case, when p = 3. Fornally, ve 


Een Sobra tie NOSE), 120}, Tori > 1,2, «se e, Ky, are 
inde- 
i k pendent renewal 
processes and N (t)= 1 (t), then u Do 
i=1 ı is 
called pooled process which is total X force casualt occus 
c: 1 rcusedbenavirour by k force 


soldiers. 


B. SCME GENERAL PROPERTIES 


1. Normal Approximation of Pooled Output 


From casualty  occurence processes in the previous 


Chapter, we assumed that there exist a mean casualty occure- 


nce time. So let u, and 4 be the nean casualty occurence 

time and corresponding ES Date i-th Y force soldier 
Beer up) ana ip) bé the mean casualty occurence tine 
Meare GOLLesponding variance of the total casualty occurence by 
k of Y fcrce soldiers. Taen the question is how are these two 


values reiated each otner ? 


From previous study, we know tnat 1) N (t) is aprroxina- 
i 


a 

tely normal with nean t/u, and variance, 6 t/u . Therefore 

i a i 
we can say that N (t) is also approximately normal with mean 
k , ] 2 3 
EN ur and variance >: G- t/yu . 2) Also assume now that 
i-1 i-1 2 a: 
either the pooled process is a renewal process or at least 


te analogous result holds, then approximately NJ (t) is norm- 


2 3 
en mean t/u(p) and variance (F (p)t/ u (p)- 


So, from 1) and 2), we must have 


t K 1 


u (p) i=l. u. (4. 1) 
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then u(p) = e ) which is the harmonic mean or the u , 
1= u 3: 
p 
2 C ES 
i 
E E E, ------ t 
E i= 1 
u (p) u, 
i 
2 2 
A 2 3 k Ci E -3 k Gi To 
en 5 =u -ER go oo (4. 
pU Gl i23 7773 1a ^ íi o3 U^) 
u, i u. 
i i 
SONT we assume that u - u and j- - TOR NAAA 222 K 
i 


which means all casualty occurences are independent and identi- 
cal, then mean casualty occurence time and corresponding vari- 


ance of the total casualty behaviour by k of Y force are as 
follows; l 


1 SA u 
ME) - ( k ) > 
2 2 
2 1 23 9 75 a EP. 
ND - (L5) AE n n) (4. 3) 
u X 


Also total number of expected casualties, E[N (t)], can 


te obtained here by fooling individual casualty occurences 


as followinrg; 


k 
E[N (t) ] * 9—, E(N (t)] = k ELM, (t)] = k M(t) (4. 4) 
p ı=1 i l 


where M(t) is expected number of casualty fron 


VMs meen, Wach is t/u tor large t. 
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let's look at an example that deals with pooling o£ 


individual renewal processes. 
Example 4.1 


Suppose of Y force artillery shells fron several 
sources are bombording the same X target. Each source hurls 
Shells at the target at a rate of 40 saells per Lour. Assume 
that the interarrival times of sheils from each single 
Source are uniformly distributed over an interval (0, a]. 
Then the question is 1) What is the probability that less 
than 825 shells will be hurled at the target from a single 
source during a 20 hour period and 2) If there are 5 sources 
from which such attacxs are launched, what is the prob- 
ability that more than 4125 ( which is 5 times 825 ) shells 
are nurled at the target in a 20 hour period ? 

To answer the first question, let \ = 40 Shells/hour 
= 2/3 shells/nin. Ihen mn 3/2 nin, when we 


assume pcisson process. But interarrival times are unifornly 


distributed over an interval (0,a]. So: Ws E X] = (0% 
a)/2= a/2, where a = 3 minutes. Then Vari X] = (a - NETA - 
min. . Now let N(t) = number oí shells hurleá in (0, t] 


from a single source. Since N(t) is normal with mean t/u and 
variance t5/w from e aeron O Dor t = 20 hours = 1200 


ninutes, the probability turns out to be as following; 


1200 
aan W / 825 - 800 ` 
au m ] 
P (N(1200)«8325] = a o | == ISS > 2 2 ee | 
\ , 1200 = 3/4 / sn 16. 33 
EY 


=:5(1.53) = 0.9370 
and for second question, 


EU (t) = number of shells hurled in (0 , tj] by i-th source 
I 


Buleresı 21,2, se, 5 .-and 


5 
N (t) E Ue) = cotal Lumber of Shells hurled in (0, 
p i= L 
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Dress sources. 
Then {N (t), t20} is the pooled process for which mean 
EPterarrival time is, 
372 B 


u 

u ze — c —-—- jrauntes. 
a k 5 10 

and corresponding variance is 


2 NS 5 ge ay 
ape) = (Sy/K ) = -;;; an 


2 
Since N (t) is normal with mean t/u(p) and WE Ata Ge E AN Up) 


for t= 1200 ninutes, 


1200 
om 3710 \ 125 
P {N (1200)<4125} = «E | -====-------- | -el mS \ 
p —\ 1200 3/100 | L 
\ 3 
EZ) 


= A 2.398897 
0.0003 


So, P(N. (1200) 24125] 


2. Poisson Process or Pooled Output 





Now let's assume that individual casualty processes are 
Ensbabilistically identical witn common underlying C.D.F. of 
RBuuadty occurence time F(z). Let F (x) be the C.D.F of 
casualty occurence time of pooled process. Then the juestion 
ES Ow iS F (x) related to F(x) ? This question can be 


answered er usas the idea of excess lire. Let Y (t) be the 
i 


excess life of an X force soldier ana Y (t) be excess life 
of total X force. Then we know that Y (t) = min( Y TE Do 

E e Aiz Y (t 7 c Fig Seles 
Then 


4 1 





n. ams e PUES res Leo aci) 


a i ne aes CCAA po 
| PW A | 
— an Jan ann e an an les ann yY onn nn nann (oo on | LM um ume Bd c due SEM Yow we ee ne eee == 
t t t a * | 
d t 
EE ! Lu | A 
er ro y ano lanar al rr yo | o - o um oz Y == == Un 
4 0 ' . 9 1 7 t 
t 0 0 t t y (EN 1 0 | 
t t q t $ | 3 ! 0 ' | 
t 0 q ' e | q 0 || 
| | | 
| MA wa: = ko so = X= = X-= === kn in are ne KRE decem ne ips 
| E Y (t) 


n 


BREMEN e ___ 


Figure 4.2 Superposition of Excess Life. 


k 
P[(Y (t)X?y) = P(Y. (t)XM (4. 5) 
p l 


Bon previous study, we know that P.D.F of E.D. for individual 
renewal processes is (1/u)í 1 - F(x)] and then we can say that 


1 
HC = F (x) ] is P.D.F o£ E.D. for the pooled casualty 
u (P p 


process. Thus, 


t 
po) 


: 
P{Y(t)>y} = | - [1 - F(x) Jdx then, 
you 
BY (t)»y) = V ---- [1 - F (x) Jäx 
p y u(p) P 


ANA Lò k 
m \ MOT EE ai rrom equation 4.5 
o oy 


If we difrerentiate with respect to 7 
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1 wo k-1 1 
ERU so hte r(x) Jax} (— ~ )f 1 - F(y¥) ] 
u(p) p u y E 


Here we xnow that u(p) = k/u, So u's are ali canceliecs out. 


Then, 


C.D.F. of pooled output is 


HB ie tn k- 1 
„mn = 1-[ m | zer) [1 = F(y) ] (4. 6) 


For example, think of the exponential distribution of 


X 
casualty occurence. Then F(x) = 1 - e A 


Then C.D.F of pooled casuaity occurerce is 


LN: k-1 -\y 
By) = 1- f \ [e ]d x ] & 
P 1 


-Ay k-1 -\ 
1- (e a < : 


BAKY 
1- e ISO EYE 


u 


which follows that pcoled output is a Poisson process with 


MENN Lere"k is tne nunber 0x Y Zorce combatants. 


Here, probably the most interesting properties of 
the pooled output refer to the "local* behaviour where xk is 
ice. Khintchine [ Ref. 13] has proposed that in the limit 
the numbers of renewals in non-overlapping intervals follow 
independent Poisson distributions, thus showing that in the 
limit the fooled output is a Poisson process. Also his 
proof does not reguire each interarrival time to be identi- 
Samy distributed. In attrition analysis in combat, it 
appears to be the same way, which means that if a large 
number of independent (and identical) casualty occurences 
are pooled, then the total force casualty occurence is 


approximately a Poisson process [Ref. 9] (Ref. 13]. Thus, 
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we can say that pooled casualty occurence by total Y conba- 
tants can be approximated by Poisson process if xe assume a 
large number of Y forces. 

"hus Einaliy, from the basic theorem or Poisson 


process, we know that if {N(t),t20} 1S Poisson process with 
rate >0, then the interarrival tines {X ,i2èl) are i.i.d. 
exponential random variable with rate A . Thus, fron 
previous results, total X casualty occurences is a Poisson 
process with rate k^, where k is the number of Y force 
te times of X casualties can be 


xponentially distributed with parameter k 
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V. COMPARISON OF RENEHAL PROCESS WITH COMPUTER SIMULATION 

In this chapter we will look at the results of  scme 
Amical combat attrition process ( Y firers with different 
rumber of passive X targets) obtained through the theoret- 
ical renewal process model which is developed in this thesis 
and computer simulaticn using random numper generator. 

As we mentioned earlier, total number or expected casu- 
aities by n Y firers turned out to be n M(t) for large t 
from the renewal theory, where we assumed that all casualty 
occurences were independent and identical. But it is clear 
that these expected casualties can not be bigger than the 
nunber of X targets. Now we will coıpare this result with 
Simulated output to see whether we can {iind any interesting 
facts such as how many X targets will be most appropriate to 
apply renewal theory to compat attrition ? or what force 
puo will be most applicable to renewal process ?, etc. 

In order to facilitate the sizulations, the assumptions 
were made that n firers shoot independently at 2 targets, 
and that the times to kill tae targets are uniformly 
distributed between [0,1]. UnenmEaoiIget 215 Killed, all 
firers shooting at it are assumed to select anew target at 
mandon from those still surviving. 

The procedure of the computer simulation is as follows 
(see Appendix A); 

1. We generate the initial kill time for each firer by 
using pseudo random nurber generator, where we assume 
that every firer's kill time is unıforıly distrikuted 
between [0,1]. 

2. Each target is indexed by integer number. 

3. We need to decide which target does the firer snoot 


at ? The target selection rule is such that first we 
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generate a random uniform nunber between [0,1], then 
we take the integer value of (random number times the 
number of targets ) plus one, which designates the 
index of target that is assigned to firer 1. Same 
procedures are repeated until all the firers are 
assigned a target, with sone targets possibly hcvina 
aore than one firer. Thus at first time step random 
number for selecting the target is generated n 
times, because we have to assign a target to eacn 
firer. These targets are killed eventually as time 
proceeds. 
Now, the target that will be killed first is the target 
which is assigned to the firer who has the minimum time to 
kill( = Tnin) among the firers, where this time to kill was 
already generated. So we now know whica target is killed. 

After the first target is killed, the number of targets 
killed is registered as one and the first time step ends. | 
Now, in the second tize step, tne new kill times between 
[0,1] are generated again for ali firers who were shooting 
Acne Zirst target. The Zirers who were still in the 
process of killing a target at tae end of the first tine 
step continue to engage the same target in the second tine 
step, where time to kill in tre second step is old kilil tiae 
( generated value at first time step ) less Tmin. 

Cnce the second target is killed, the second time step 
ends and the number of targets kilied is now registered as 
two. Now we proceed to the third time step, "here same 
target selection rule and firing rule are applied until all 
the targets are killed. 

After all the targets are killed, each time step length 
( which is the same as the time to kill the target in each 
tine ster ) is stored, and this is one replication (see 
Appendix E). To obtain tne expected value, 40 replications 


rhe simulation were conducted and o data sets 
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(n3, m=6), (n=3,m=10), (c=3, m=14), and (n=10, m=6), (n=10, 
m=10), (n=10, m=14), were examined. Arter 40 replicaticns , 
the data were evaluated to find the expected number of kiils 
Besough order statistics. This expected number of kills was 
plotted against tine. The results for n=3 and n=10 are 
ENnown in Fig 5.1. 

Ihe theoretical results can be obtained through two 
methods, whith are discrete approximation and simulation. 
In this thesis, we used discrete approximation based on the 
renewal function modelled at equation 3.7. It's formula is 


as follows; 


Mer 


M(t) = 51 1 € H(t-j) ] P. 
J 


j=1 
where P represents tne probability distribution 
O Ene 


Thus, expected  nuaber of kills obtained from the reneval 
process model using discrete approximation (see appendix C 
EBorsprogram and data, where left colunn refers to the kill 
tine and right coluna refers to the expected number of kills) 
is adse Flottea against time in fig 5.1. 

It was round through both cases that the expected number 
of kills ty the renewal process model may only pe justified 
Exrougb the assunption of an infinite number of targets 
because as m gets bigger, the expected kiil value approacnes 
tne theoretical result. Also it was discovered tnat at the 
beginning of contrat, the attrition process almost followed 


the theoretical renewal process modei. 


47 





14 
wa 2° M= 10 
6 


3 FIRERS 
Mz 
2 
TIME 
10 FIRERS 
14 
M - 10 
2 


THEORETICAL VS SIMULATED RESULTS 
M= 


MOT  ~0E T 9 of 0% m 
SIA JO W3GN^N Q3193dX3 S115 JO YIBNNN 03193d4X3 


TME 


Expected Kills From Theoretical And Simulated. 


Figure 5.1 


48 





EX SUMMARY 


Since Larchester proposed his model, many studies have 
tried to model the ccmbat attrition process mathematically. 
Recently, the main focus has been on calculatior of numer- 
ical values for Lanchester attrition-rate coefficient. This 
numerical determination stands at tne heart of casualty 
assessment in such models. There are three reasons for the 
importance of this attrition-rate  coefricient methodology. 
First,  Lanchester-type models are used in various U.S.Army 
and D.O.D. planning activities nore widely than ever before. 
Second, Lanchester attrition-rate coefficient is a tasic 
element of any  Lanchester-type combat model, and that 
attrition-rate coefficient reflects tnae effective applica- 
tion of firepower. Finally, significant new developments 
have occured in nethcdology for developing more tactically 
realistic Lanchester attrition-rate coefficiert and these 
important results have not been accessible to a very wide 
audience. In particular, a new approach for developing more 
realistic and mathematical Lanchester attrition-rate coeffi- 
cient, that of computing combat attrition by using a renewal 
process, was proposed py Taylor recently. This thesis began 
with the assumption that the combat attrition process was a 
Kind of renewal process, and that tines between casualties 
Mememaei.d. Thus, considering the problem from the stand- 
point cê renewal theory, many new random variables were 
included and various distribution functions studied in order 
to understand better this type of casualty occurence 


paradign. 
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This thesis only dealt with tne case of homogeneous 
forces which were using the same acquisition mode to attrit 
the targets, which were passive X forces. Thus we Could use 
an ordinary renewal process to study the behaviour oí casu- 
alty occurences. We emphasized the probabilistic function 
of each r.v., so once we found the distribution of interar- 
rival times of casualties, casualty analysis could be Gone 
ty plugging distribution functions into the equation devel- 
oped. HER ERE AIlSEE1TBUution of interarrival times of 
casualties was the critical factor in continuing the casu- 
alty analysis- 

The final conclusion is that total casualty occurence in 


combat can be assumed to be a Poisson process, ard thus 


ip 


times between casualties are exponentially distributed if we 
assume a large number of combatants. 

We have dealt here only with a passive X force( i.e. one 
that does not shoot rack ), perhaps an unrealistic assuanp- 
Hon. But still tnis study permits Letter understanding or 
Senbat attrition, through the use of probability theory, 


than was previously available. 


B. SUGGESTIONS FOR FURTHER STUDY 


In this thesis, we have only dealt with passive 
targets. A rore complex formulation may be rezuired for real 
compat, since enemy targets shoot back and cause firer 
ERSETEI tion. 

In addition, a delayed renewal process could be applied 
for the case where different acquisition modes are used for 
the first shooting and for each subsequent shootings. À 
reinforcement of the rorces Could also bo Included: 
Casualty analysis could then be done using the theory of 
aiternating renewal processes, once we know the rate of 


reinforcement. 
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FORTRAN PROGRAMMING FOR SIMULATING TIE KILL TIME (N = 3) 
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FUNCTION 


DISCRETE APPROXIMATION FO RENEHAL 
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